Crossover in the scaling of island size and capture zone distributions 
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Simulations of irreversible growth of extended (fractal and square) islands with critical island 
sizes i = 1 and 2 are performed in broad ranges of coverage 6 and diffusion-to-deposition ratios R 
in order to investigate scaling of island size and capture zone area distributions (ISD, CZD). Large 
9 and small R lead to a crossover from the CZD predicted by the theory of Pimpinelli and Einstein 
(PE), with Gaussian right tail, to CZD with simple exponential decays. The corresponding ISD 
also cross over from Gaussian or faster decays to simple exponential ones. For fractal islands, these 
features are explained by changes in the island growth kinetics, from a competition for capture of 
diffusing adatoms (PE scaling) to aggregation of adatoms with effectively irrelevant diffusion, which 
is characteristic of random sequential adsorption (RSA) without surface diffusion. This interpreta- 
tion is confirmed by studying the crossover with similar CZ areas (of order 100 sites) in a model 
with freezing of diffusing adatoms that corresponds to i = 0. For square islands, deviations from 
PE predictions appear for coverages near 6 = 0.2 and are mainly related to island coalescence. Our 
results show that the range of applicability of the PE theory is narrow, thus observing the predicted 
Gaussian tail of CZD may be difficult in real systems. 

PACS numbers: 68.43.Hn, 05.40.-a, 68.35.Fx , 81.15.Aa 



I. INTRODUCTION 

The theoretical study of thin film and multilayer 
growth attracted much interest in the last decades mo- 
tivated by the increasing number of experimental tech- 
niques and applications^— . Film morphology is strongly 
connected to the initial stages of its growth, where is- 
lands of various shapes may be formed. The submono- 
layer regime, in which a single incomplete layer is being 
formed, was modeled by various authors and is still the 
focus of much debatei^— . The theoretical approaches 
usually try to explain the island size distributions (ISD) 
and the capture zone distributions (CZD), with a cap- 
ture zone (CZ) defined as the area in which a diffusing 
adatom is more likely to attach to a given island than to 
any other one. The models were originally proposed for 
atomic epitaxy and recently extended to growth of other 
materials, such as organic molecule islands, colloidal epi- 
taxy, and graphene epitaxy^—. 

For irreversible island growth, the most recent advance 
in the field is the theory of Pimpinelli and Einstein (PE) 8 , 
which proposed that CZD are described by the Wigner 
surmise (WS) from random matrix theory^. The basis 
of PE theory is the competition of CZs for aggregation 
of diffusing adatoms. Deviations in the peaks of simu- 
lated CZD and the WS were reported for point islands 
on surfaces^. However, a recent work showed excellent 
collapse of most CZD with the WS after suitable rescal- 
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ing, for point, fractal, and square islands, and showed 
that all CZD had the universal Gaussian decay predicted 
by PEi 2 -. Indeed, that theory obtains the WS from a 
Langevin equation, which is expected to apply only for 
large CZs. Moreover, the tail of the ISD was predicted 
by a scaling approach that assumes the Gaussian decay 
of CZD and accounts for the island shaped 2 -. 

On the other hand, failure of the predicted decay of 
CZD was observed for point islands in one-dimensional 
lattices, where a decay as exp(— x 3 ) is found 2 ^— . This 
decay is also predicted analytically in two dimensions 
under a circular-cell approximation for CZs 2 ^. For the 
more realistic cases of extended islands, there are ad- 
ditional reasons to expect limitations of the PE theory. 
First, as increases, the size of island-free regions de- 
creases, enhancing the correlations between stable islands 
and consequently changing their competition dynamics. 
Moreover, if the diffusion-to-deposition ratio R = D/F is 
small, the diffusion lengths of the free adatoms are small, 
facilitating formation of new islands instead of capture by 
the existing ones. 

This work is devoted to study the crossover from PE 
submonolayer growth (Gaussian CZD and related forms 
of ISD) in those conditions. 

Simulations of irreversible growth of fractal and square 
islands with critical island nucleus of sizes i = 1 and 2 
are performed with broad ranges of 9 and R. There is 
a variety of shapes of CZD and ISD for small CZ areas 
and small islands (left tails of CZD and ISD). On the 
other hand, the right tails of CZD and ISD always show 
a crossover to simple exponential decay for small R or 
large 0. In some cases, it restricts the PE prediction to 
a narrow ranges of those variables. 
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Those deviations are related to the decrease of the size 
of free regions between the islands, which become very 
close and strongly correlated, ruling out the (mean-field) 
PE theory that neglects island correlations. For small R, 
the exponential tail is an effect of small diffusion lengths 
of adatoms before aggregation, typical of random distri- 
butions of non-mobile adatoms, which is the random se- 
quential adsorption (RSA) proble m 24 i 25 . As 9 increases, 
for any R, fractal islands become large, thus most of the 
atoms are deposited near their branches and rapidly ag- 
gregated. Again, adatom diffusion lengths are small and 
RSA scaling appears. A CZ area reduction to less than 
100 sites is characteristic of this crossover region. This 
is confirmed by results of a model with i = where dif- 
fusing adatoms may stop moving without aggregating to 
an existing island (adatom freezing) 6 . The crossover in 
square island models as 9 increases and the islands begin 
to coalesce is interpreted along the same lines. 

The study of ISD and CZD is also important for ex- 
perimental work 2 ^— , particularly to determine the crit- 
ical island size that reveals basic features of the aggrega- 
tion processes. At this point, we note that recent works 
compared CZD and ISD with fitting curves including 
the WS, such as para-sexiphenyl island growth on dif- 
ferent substrate a 30 ' 31 , Cu deposition with impurities 32 , 
pentacene island growth with impurities^ 7 -, InAs quan- 
tum dot growth on GaAa^, and Ceo deposition on SiC-2 
films^. 

The rest of this work is organized as follows. In Sec. 
ITU we present the growth models and summarize previ- 
ous theoretical results. In Sec. IIII| we discuss the scaling 
of CZD and ISD of fractal islands for broad ranges of 
coverage and diffusion-to-deposition ratio, with special 
emphasis on deviations from PE predictions. In Sec. IIV( 
the discussion is extended to square islands models. In 
Sec. [V] we analyze a model where diffusing adatoms 
may freeze, which presents a crossover from regimes with 
relevant and irrelevant adatom diffusion. Sec. IVII sum- 
marizes our results and presents our conclusions. 



II. MODEL DEFINITION AND THEORETICAL 
APPROACHES 

We consider models of atom deposition and diffusion 
on a square lattice with irreversible aggregation to is- 
lands larger than a critical size i. Atoms incide at ran- 
domly chosen lattice sites with rate F (number of de- 
posited atoms per unit time per lattice site). An atom 
adsorbs at the incidence site only if it is empty, thus 
growth is restricted to a single layer. An adatom diffuses 
with coefficient D (number of random steps to neighbor- 
ing sites per unit time) if it is not aggregated to a stable 
island. The critical nucleus i is the maximum size of a 
non-stable island, i. e. an island from which atoms can 
detach. Islands of size i + 1 or larger are stable. 

The diffusion-to-deposition ratio is defined as R = 
D/F. The coverage 9 is defined as the ratio between 



the number of adsorbed atoms and the number of lattice 
sites. 

In fractal island models, each atom permanently ag- 
gregates at the site where it collides with a stable is- 
land. Since there is no relaxation of the aggregated 
atoms, branched shapes are produced, resembling those 
of diffusion-limited aggregation— 

In square island models, after aggregation to a sta- 
ble island, the atom instantaneously relaxes to a position 
that preserves the compact (square) shaped. However, 
when two neighboring islands coallesce, they evolve in- 
dependently: an atom aggregated at the border of one 
of the coallesced islands relaxes to a site in the border of 
that island to preserve its square shape. 

In the case of point island models, all atoms of an 
island aggregate on a single lattice site. However, this 
case is not studied in this work, since our aim is to study 
features of extended islands. 

It is reasonable to assume that the probability density 
of CZ area x follows the scaling form 
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where / is a scaling function. An equivalent scaling form 

applies to the density of islands of size s, Q(s). A recent 

work 12 proposed the scaling of CZD with the variance 

2 1/2 

a x — (x — (x)) as 



Of N 1 ( X ~ ( X ) 

P(x) = —g[ 
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This procedure is inspired in the successful application to 
roughness distributions of thin film growth model o 36 ^ 37 . 
An equivalent form also applies to Q(s). 

In 1995, Amar and Family (AF) proposed the empirical 
formula for the ISD 6 



fi (u) = CiU 1 exp (~ia l u 1/ai ^j , 



(3) 



where u = s/(s) and Ci and are normalization con- 
stants. This formula is widely used to fit experimental 
data— r— . Subsequently, Mulheran and Blackman (MB) 
proposed to relate the ISD to distributions of areas of 
Voronoi polygons and provided results close to those of 
point island models 7 . In all those works, comparison with 
simulation or experimental data focused on the peaks of 
the distributions. 

A recent work by Korner et al 3 ^ proposed mean-field 
rate equations accounting for coverage-dependent cap- 
ture numbers and obtained good agreement with the 
peaks and the tails of the ISD for fractal islands. How- 
ever, the theoretical ISD were obtained by numerical in- 
tegration of those equations, thus the tail decays were 
not analyzed. 

Recently, a significant advance was the proposal^ that 
CZD are described by the WS 
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where z = x/ (x), ft = § (i + 1), d is the substrate dimen- 
sion (d — 2 in the present work) and the parameters ap 
and bp are determined by normalization conditions. The 
PE theory is justified by a phenomenological argument 
that the CZD can be extracted from a Langevin equa- 
tion representing the competition of neighboring islands 
for adatom aggregation. 

Ref^ showed that the CZD scaling using Eq. (O 
agrees with the WS of Eq. Q for large values of R (typi- 
cally R > 10 7 ) and small coverages, in point, fractal, and 
square island models. In all cases, the Gaussian right 
tail predicted by PE [Eq. (Q}] was present. For point 
and square islands, some deviations in the left tails and 
in the peaks of the CZD appeared. However, deviations 
for small islands and small CZs are expected because the 
continuous approach of PE applies to large island and 
large CZ kinetics. 

Refi 2 - considered /3 = i + 1 in the WS, as originally 
suggested 8 . The successful comparison is probably due to 
the adopted rescaling and contrasts to the proposal (3 — 
i + 2 of Refsi^r— . The Gaussian decay also differs from 
the exp[— x 3 ) decay predicted in Refi^2 for point islands, 
probably due to the use of a circular-cell approximation 
for CZs in that work. 

A scaling approach predicts the decay of the right tails 
of ISD from the Gaussian tail of the CZEU^. For fractal 
islands, the decay is exp (— s 4 / Df ) « exp (— s 2 36 ), where 
D F ss 1.694 39 is the fractal dimension of DLA clusters; 
for square islands, the ISD decay is also Gaussian. These 
results also agree with simulation data for large R and 
small coverages 1 ^. 

PE theory applies to systems where distant and large 
islands compete for aggregation of diffusing adatoms. 
The opposite situation is that of totally non-mobile 
adatoms, which is the random sequential adsorption 
(RSA) problem without diffusio n 24 ' 25 . In this case, atoms 
irreversibly stick to the site where they are adsorbed. ISD 
and CZD have simple exponential decays 4 ^ as 

Q(«)~ exp (-«/(*»• (5) 

III. DISTRIBUTIONS FOR FRACTAL ISLANDS 
A. Simulation results 

We performed simulations in lattices of very large lat- 
eral size, typically L = 2048, with R ranging between 
10 4 and 10 9 , coverages up to 9 = 0.4, and critical nu- 
cleus i = 1 and i = 2. For i = 2, no energy barrier is 
considered for the dissociation of islands with two atoms. 

In Fig. la, we show scaled CZD for R = 10 9 and 
various coverages, to be compared with the WS. The data 
for small coverage (9 = 0.1 in Fig. la) is well fit by the 
WS, as previously shown in Refji 2 . However, deviations 
clearly appear for the largest coverages (9 = 0.3 and 0.4) 
and a simple exponential decay of the right tail is found, 
in contrast to the Gaussian decay of the WS. 
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FIG. 1: (Color online) Scaled CZD (a) and ISD (b) for fractal 
islands with i = 1, R = 10 9 and several coverages. The solid 
line is the WS with ft = 2. Dashed lines in all plots are guides 
to the eye. 



In Fig. lb, we show ISD for the same set of parameters. 
For small coverages, the right tail decay is slightly faster 
than a Gaussian, as discussed in Sec. |H](see also Refi 1 ^). 
However, it also crosses over to simple exponential decay 
for the largest coverages (9 = 0.3 and 0.4). 

The same crossover is found as R decreases, for fixed 
9. In Fig. 2a, we show CZD for 9 = 0.10 and several 
values of R. The CZD for R = 10 8 agrees very well with 
the WS. However, the CZD for smaller values of R devi- 
ate from the WS and these deviations are enhanced for 
decreasing R. The simple exponential decay is found for 
large CZs for R < 10 5 . Despite the deviation in the tails, 
for R = 10 6 , the peak of the CZD still shows reasonable 
agreement with the WS (inset of Fig. 2a). For smaller R, 
the peaks also deviate from the WS, compensating the 
large difference in the tails. 

In Fig. 2b, we show CZD for 9 = 0.3 and several values 
of R. For R = 10 9 , the CZD agrees with the WS in two 
orders of magnitude around the peak [a z P (z) ~ 10~ 2 ' 5 
to 10~°' 5 ]. For smaller values of R, the fit covers a smaller 
region of the scaled CZD. The simple exponential decay 
of the right tail is observed for large CZs in all cases. 
The corresponding ISD are shown in Fig. 2c and also 
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FIG. 2: (Color online) (a) Scaled CZD for fractal islands with 
i — 1, 9 — 0.10 and several ratios i?. The inset in (a) shows 
the same data of the main plot in linear-linear scale. Scaled 
CZD (b) and ISD (c) for fractal islands with i = 1, = 0.30 
and several ratios R. The solid line is the WS with /3 = 2. 
Dashed lines in all plots are guides to the eye. 



have that decay. 

Note that coverage 9 = 0.3 is typical of many exper- 
imental works because islands become sufficiently large 
and facilitate microscopy investigation with a focus on 
island statistics. Moreover, R < 10 5 is typical of low 
temperatures, where island relaxation is difficult and the 



fractal island model is realistic. These features suggest 
that experimental observation of PE predictions [partic- 
ularly the Gaussian tail of CZD and the related form 
of the ISD (Sec. fll])] will be difficult when large fractal 
islands are grown. 

Fig. 3 shows images of the submonolayers grown with 
R = 10 7 and R = 10 9 . For 6 = 0.1 and R = 10 9 , 
the clusters are distant from their neighbors, i. e. there 
are large areas between them. The PE theory actually 
applies to this situation^. However, for R = 10 7 and the 
same coverage, islands are much smaller, the distance 
between them decreases and the free areas around them 
are smaller. Under these conditions, the hypothesis of 
PE theory begin to fail. For 9 = 0.1 and R = 10 6 , 
these effects are enhanced, which explains deviations of 
the CZD from the WS, as shown in Fig. 2a. 



B. Interpretation of results 

First we consider the deviations from PE scaling ob- 
served for small values of R. 

In this case, the average adatom diffusion length (Id) is 
very small, as well as its average diffusing time (to)- Is- 
land growth mainly proceeds by incidence of new atoms 
very close to existing islands or very close to a diffus- 
ing adatom. If the lattice size is rescaled by (Id) and 
time is rescaled by (tn), then the islands of the rescaled 
lattice will grow by attachment of atoms deposited at 
neighboring sites, i. e. non-diffusion atoms. This leads 
to distributions with right tails (large islands) similar to 
RSA (Eq. El). 

Another important aspect to explain the deviations 
from PE scaling for small R is the small average island 
size and small average CZ area, since that theory follows 
from continuous approaches suitable for large CZs. 

This interpretation helps to quantify the deviation of 
the DCZ from the WS as R decreases. 

The deviation occurs for a CZ area x* in which the 
WS decay exp —(x/(x)) 2 (Eqs. [T]and|3]), matches the 

simple exponential decay of Eq. [5j Thus x* ~ (x) . In 
the crossover, islands are large and very close (see e. g. 
Fig. 3), branched but not fractal, thus the CZ area x* 
is of the same order of the average island size. The de- 
pendence of this size on R and 9 in the steady state of 
island growth is predicted by rate equation theorji^ir— 

as (s) - iT/(i+2)0-(i+l)/(i+2)_ Thig giveg 



R i/(i+2)Q-(i+l)/(i+2)_ 



(6) 



The crossover area was determined in the simulated 
CZD as the point in which the tail changes its concativity 
in a log-linear plot. Error bars were estimated as the 
whole range of x that covers the crossover region. Fig. 
4 shows the scaling of x* with R, for fixed 9 = 0.3, with 
linear fits giving scaling exponents \ = 0.33 ± 0.02 for 
i = 1 and \ = 0.47 ± 0.03 for i = 2. This result is 
in good agreement with the predictions x* ~ i? 1 / 3 and 



5 







, ^ _ I. 




^^^^ 






FIG. 3: (Color online) Fractal islands for several coverages with R 
size 200 sites are cut from a system of size L = 2048. 



10 7 (top) and R = 10 9 (bottom). The panels of lateral 




FIG. 4: (Color online) CZ crossover area x* as a function of 
R, for 8 = 0.30 and critical islands sizes i = 1 (squares, red) 
and i = 2 (triangle, blue). The error bars are of the same 
order of the point size. Dashed lines are least squares fits of 
the data. 



x* ~ R 1 / 2 of Eq. © for i = 1 and 2, respectively. The 
scaling of the average island size at the crossover, s* , 
is analogous to a;*, and was also confirmed numerically 
(results not shown). 

Another interesting conclusion from this quantitative 
analysis is that the rate equation theory is still a reliable 
approximation for average quantities, despite the changes 
in the main kinetic mechanisms determining the evolu- 
tion of large CZs and large islands (responsible for the 
crossover in the right tails of ISD and CZD). 

Now we consider the deviations from PE behavior with 
increasing coverage. For instance, Fig. 3 shows that 



the islands are very close (so, strongly correlated) for 
R = 10 7 and 9 = 0.2. They grow mainly by random 
aggregation of mass deposited near its branches, and not 
by capturing diffusing atoms that incided in island-free 
areas. Thus, there is significant increase in the widths 
of island branches when 9 evolves to 0.3 and 0.4, instead 
of increase of the length of the branches (in other words, 
islands become less ramified). 

This situation may also be viewed as a case of small 
average adatom diffusion lengths because the free area 
available is small. Again, rescaling of the lattice by (Id) 
and of the time by (to) indicates that most atoms will 
aggregate immediately after incidence. This corresponds 
to RSA scaling. 

It is worth to recall that an attempt to deposit an 
atom on an already occupied site is rejected in our mod- 
els, and a new site is randomly chosen. However, in many 
simulation works, that atom is randomly aggregated in 
the island periphery. For low coverages, this condition 
is irrelevant, but for large coverages the second rule in- 
creases the number of static adatoms aggregated to is- 
lands. Thus, the crossover effects would be enhanced if 
this (more realistic) mechanism is used. 



IV. DISTRIBUTIONS FOR SQUARE ISLANDS 

We performed simulations of the square island model 
in the same range of parameters of the fractal island sim- 
ulations. 

In Figs. 5a and 5b, we show CZD for R = 10 9 and 
various coverages, in logarithmic and linear plots, respec- 
tively. The data for small coverage (9 = 0.1) is well fit- 
ted by the WS, but deviations are large for 9 = 0.25 
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be noted in the distributions of Figs. 5a-c. 

The first one is the presence of many small islands for 
all coverages, which was already noted in Ref. 44 . For 
9 = 0.4, the density of small islands (s — 2,3,4) is 
comparable to the density of islands of the average size 
(s fx 1100). This certainly affects all results from contin- 
uous approaches, which is the case of the PE theory, and 
justifyes the deviations between the CZD and the WS for 
relatively small coverages. 

Another difference is the formation of a shoulder in the 
CZD as 9 increases (Fig. 5b). This is also related to the 
presence of small islands. Despite their small size, their 
CZ areas are not small; instead, they are only slightly 
smaller than the most probable CZ areas (peaks of CZD). 

Results for fixed coverage and decreasing R show the 
same crossover of the fractal islands (Figs. 2a and 2b). 
However, the simple exponential decay of the right tails 
(of CZD and ISD) appears for larger values of R, i. c. 
the Gaussian CZD is observed in a narrower range for 
square islands. 

Fig. 6 shows images of the submonolayers grown with 
R = 10 7 and R = 10 9 . For any coverage, the free areas 
around the square islands are larger than those around 
the fractal islands. Thus, one should expect PE theory to 
apply for larger coverages with square islands. However, 
the opposite occurs: deviations of the CZD from the WS 
appear for smaller coverages when compared to fractal 
islands. 

The formation of small islands explains the deviations 
in the left tail and in the peak of CZD. On the other 
hand, the change in the right tail scaling (from Gaussian 
to simple exponential) is not explained by this feature, 
but is related to island coalescence. For 9 = 0.2, only 
two-island coalescence can be seen in Fig. 6 and is sig- 
nificant for R = 10 7 . For 9 — 0.4, coalescence of many 
islands is noticeable for both values of R. Now suppose 
that the system is rescaled by a factor (s), so that an 
average island of the original system is converted to a 
single occupied site. After rescaling, the largest islands 
are formed by coalescence of some randomly distributed 
sites, which were the coalesced islands in the original sys- 
tem. Thus, ISD and CZD have right tails as in random 
distribution of occupied sites, i. e. RSA. 



FIG. 5: (Color online) Scaled CZD (a,b) and ISD (c) for 
square islands with i = 1, R = 10 9 and several coverages. 
The solid line in (a) and (b) is the WS with j3 = 2. Dashed 
line in (a), (c) are guides to the eye. 



and larger coverages. Again, the Gaussian decay of the 
right tail of the CZD (PE scaling) crosses over to simple 
exponential decay as 9 increases. 

In Fig. 5c, we show ISD for the same set of parameters, 
showing the crossover from the Gaussian decay (small 
coverage) to the simple exponential (large coverage) . The 
inset of Fig. 5c highlights the small island densities. 

Important differences from the fractal island case can 



V. MODEL WITH DIFFUSION AND FREEZING 
(i = CASE) 

This model was introduced by Amar and Family^ to 
represent growth with critical island size i = incorpo- 
rating diffusion effects. That model is different from the 
random sequential adsorption (RSA), in which atoms ir- 
reversibly stick to the site where they are adsorbed, with- 
out diffusion^ 4 ^. 

In the model, the adatoms diffuse with coefficient D 
and stop moving (freeze and nucleate a new stable is- 
land) with rate R s = rD. with r < 1. The freez- 
ing rate depends on an extra activation energy Ep as 
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FIG. 6: (Color online) Square islands for several coverages with R = 10 7 (top) and R = 10 9 (bottom). The panels of lateral 
size 200 sites are cut from a system of size L — 2048. 



r = cxp (—Ep/ksT), which may represent the effect of 
surfactants or impuritie a 45 ' 46 . Adatoms also stop mov- 
ing when they become nearest neighbors of stable is- 
lands, producing islands with fractal shape. However, an 
adatom does not stop moving when it meets one or more 
diffusion adatoms, which is the main difference from the 
models with i > 0. 

We performed simulations of this model in lattices of 
lateral size L = 2048, with R between 10 9 and 10 11 , 
coverages up to — 0.3, and r ranging from 10~ 2 to 
10~ 5 . 

In Figs. 7a and 7b we show scaled CZD and ISD, 
respectively, for different freezing rates and fixed R and 
9. For r = 10~ 5 , the CZD shows a reasonable agreement 
with the PE theory for i = (WS with [3 = 1). However, 
a large deviation is observed for r = 10~ 3 , where the tail 
decay is slower than a Gaussian (see inset of Fig. 7a). 
The ISD have simple exponential decay for r = 10~ 3 and 
large s, again typical of RSA. 

We notice that, in this model, islands are produced 
only by adatoms that stop moving. Thus, the number of 
islands increases with the rate r, decreasing the CZ areas 
and the size of the free regions between the islands. 

The crossover in CZD can be explained by a compari- 
son of CZ areas with the model with i = 1 (Sec. 

The average CZ area for r = 10~ 5 and 8 = 0.1 is 
(x) « 880. Comparing with fractal islands grown with 
2 = 1 and the same coverage, we note that this value is 
between the CZ area for R = 10 7 ((x) « 500) and for 
R = 10 s ((x) k 1200). Both CZD agree with the WS 
with f3 — 2, as shown in Refi 1 ^ (see also Fig. 2a for 
R = 10 8 ). The case R = 10 7 is also illustrated in Fig. 3 
and show large areas between the islands. 

On the other hand, the average CZ area for r — 10~ 3 



and 6 1 = 0.1 is much smaller: (x) « 88. This value 
is close to the average CZ area for fractal islands with 
i = 1, 6 = 0.1, and R = 10 4 : (x) w 64 (for R = 10 5 , 
we find (x) ~ 120). The corresponding CZD is shown in 
Fig. 2a and has significant deviation from the WS with 
/3 = 2. 

The deviations from PE scaling for small R in fractal 
islands are related to the decrease of the average CZ area 
and the corresponding decrease in the available area be- 
tween the islands. The crossover occurs for an average CZ 
area near 100 lattice sites, which apparently independs 
on the value of i, reinforcing the geometric argument. 

The case r — 1 is equivalent to the RSA, since adatoms 
freeze with the same rate that they move to neighbor- 
ing sites. Consequently, our results consistently illustrate 
a crossover from a regime of island growth from diffus- 
ing adatoms to a regime of island formation by approxi- 
mately static adatoms. 



VI. CONCLUSION 

We simulated irreversible growth of fractal and square 
islands with critical island size i = 1 and 2, for several 
values of diffusion-to-deposition ratio R and coverage 9, 
and simulated a model with i = 0, with adatom diffu- 
sion and freezing. ISD and CZD were calculated and the 
crossover from PE theory predictions to different distri- 
bution shapes was analyzed, with a focus on the right 
tail decays. 

For fractal islands, small 8 and large values of R lead 
to Gaussian CZD, which is consistent with the PE the- 
ory. ISD have left tails slightly faster than Gaussian. 
However, as the island density increases (due to increas- 
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FIG. 7: (Color online) Scaled CZD (a) and ISD (b) for i = 
model with R = 10 11 , (9 = 0.10 and r = 10~ 3 (squares, red) 
and r = 10 _J (triangles, blue). The solid line is the WS with 
(3 = 1. The insets show the same distributions with abscissas 

2 in (a) and [<£^£»] " in (b). 



aggregating to a stable island, which is interpreted as a 
crossover to the problem of random distributions of static 
atoms (RSA problem). For fixed coverage, estimates of 
the island sizes and CZ areas at the crossover have R- 
scaling consistent with rate equation theory for i = 1 
and i = 2. 

For square islands, small 9 and large values of R lead to 
Gaussian CZD and ISD. Compared to the fractal island 
case, the crossover is observed for smaller values of 9 and 
larger values of R, despite the larger area between neigh- 
boring islands. In the square island case, the deviations 
from PE predictions is partially atributed to the presence 
of many small islands, whose densities are of the same or- 
der of the (much larger) average islands. However, the 
crossover from Gaussian to simple exponential decay in 
CZD and ISD is mainly related to coalescence of islands, 
since large islands and large CZs are effectively formed 
by random distribution of (static) rescaled islands. 

In the model with diffusion and freezing, island shape 
is also fractal. For small values of the freezing rate r, 
islands are separated by large areas and the CZD obeys 
the PE theory for i — 0. As r increases, the area between 
the islands become very small, so that CZD and ISD get 
a simple exponential decay, similar to the fractal islands 
with i = 1. 

The Gaussian tail of CZD and related forms of ISD 
are found only for small coverages (usually 9 < 0.3) and 
for large values of R, thus their experimental observation 
is probably difficult. Moreover, large values of R are 
characteristic of high temperatures, where the hypothesis 
of a critical (and small) island size may fail. This strongly 
suggests extensions of the theoretical and simulational 
study to reversible island growth and to systems with 
cluster diffusion, which were focus of recent work 47 . 



ing 9 or decreasing R), the area between neighboring is- 
lands becomes very small. This strong correlation be- 
tween islands leads to a crossover to simple exponential 
decays of CZD and ISD. For small R or large 9, the de- 
posited adatoms have very small diffusion lengths before 
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